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Abstract. We have performed MD simulations of a highly charged colloid in a solution of 3:1 and addi- 
tional 1:1 salt. The dependency of the colloid's inverted charge on the concentration of the additional 1:1 
salt has been studied. Most theories predict, that the inverted charge increases when the concentration 
of monovalent salt grows, up to what is called giant overcharging, while experiments and simulational 
studies observe the opposite. Our simulations agree with the experimental findings and shed light onto the 
weaknesses of the theories. 

PACS. 64.70.pv Colloids - 41.20.Cv Electrostatics; Poisson and Laplace equations, boundary- value prob- 
lems - 87.10.Tf Molecular dynamics simulation 



1 Introduction 

In this work we study the phenomenon of overcharging on 
the basis of the simple restrictive primitive model. When 
a highly charged colloid (or macroion) is put into a solu- 
tion that contains multivalent counterions, its charge can 
become overcompensated, such that the effective charge 
of the colloid-ion-complex becomes oppositely charged. In 
the literature, this phenomenon is called charge reversal, 
charge inversion or overcharging. 

Experimentally, the phenomenon has been demonstrated 
by electrophoresis experiments, where a reversed mobility 
has been observed for charged colloids in a solution that 
contained multivalent ions (for a review, see [1]). Only re- 
cently, experiments were able to directly measure charge 
reversal in a system of a solution of multivalent salt at a 
charged wall [HIS]. 

From a theoretical point of view, the standard Poisson- 
Boltzmann theory is not able to predict charge inversion. 
A number of theories such as integral equations [4,5,61 
[7], field theoretic approaches [8jl9j, and others, see for ex- 
ample the references in [TCTlfTT] . have been put forward 
to describe this effect. Most popular amongst these are 
the simple one-component plasma (OCP) theories, which 
originated in the Wigner-Crystal theory by Shklovskii p!2l 
[13] , and have some descendants by Levin [Tl] and a recent 
one by Pianegonda et al [TS]. The basic idea of these the- 
ories is, that on the surface of the colloid, the multivalent 
ions form a Wigner crystal or at least a strongly correlated 
liquid. Therefore, the charge can be overcompensated by 
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the multivalent ions. The rest of the system is considered 
in a Poisson-Boltzmann-like fashion. Although these theo- 
ries are pretty simple, they give good predictions in many 
cases. 

Simulations have confirmed the phenomenon of over- 
charging P^[r7l [TS |ll9ll20| . and they are mostly consistent 
with the OCP theories in the strong-coupling regime, which 
is close to the Wigner-crystal (T — 0) limit. A notable ex- 
ception is the work by Messina et al [2T] , which shows that 
overcharging is also possible in situations that are very far 
from the strong coupling limit. Here, the overcharging was 
attributed to mostly an entropy-driven effect. 

Realistic and in particular biological systems, very of- 
ten contain a significant amount of "normal" monovalent 
salt. Therefore it is interesting to see what happens to the 
charge reversal when 1:1 salt is added to the system. 

Experimentally, with growing concentration of mono- 
valent salt, the inverted charge seems to decrease [T]. In 
contrast to that, most OCP theories predict a monotonous 
increase of the inverted charge when more salt is added. In 
an extension to Shklovskii 's original theory, Nguyen, Gros- 
berg and Shklovskii [22j,23j even predict so-called "Giant 
overcharging" for high concentrations of monovalent 
salt, where the effective charge is larger than the bare (neg- 
ative) charge of the colloid. Levin's theory [2] predicts a 
more moderate growth, while Pianegonda, Barbosa and 
Levin |15] have a more detailed view: for low salt con- 
centrations, an increase for growing salt concentrations is 
predicted, which reverses into a decrease for higher con- 
centrations. 

Not many simulations have been performed of that 
situation, however, all of them seem to support the ex- 
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perimental finding, that the amount of reversed charge 
decreases with increasing salt content. Tanaka and Gros- 
berg have performed simulations of the electrophoresis of 
a colloid in salt solution in a regime of very high sur- 
face charge ^U\. Martin-Molina [^[^ et al simulate the 
overcharging at a planar wall in salty environment un- 
der more realistic conditions, which should be valid in the 
limit of large colloids. Diehl and Levin [26J have performed 
Monte-Carlo simulations of a spherical colloid, where the 
^-potential is measured. 

This work complements these simulations by analyzing 
the case of a spherical colloid in salty environment in the 
experimental parameter regime, and comparing the results 
to the various assumptions made in the OCP theories. 



2 System 

The model system we are considering is a single, spherical, 
highly charged colloid in a solution of asymmetric 3:1 salt 
and additional 1:1 salt in a cubic box with periodic bound- 
ary conditions and a side length of L = 225.8 A, which 
corresponds to a spherical cell with a radius Rq — 140 A. 

The simulation box contains four particle types that 
are defined via their diameter dj and their charge qi : 

- 1 charged colloid (Qq = -300 eo, do = 100 A) 

- 200 (-1-3) counterions ((7(+3) = 3eo, cJ(+3) — 4 A) 

- (300 + TT-sait) (-1) colons (g(_i) = -1 eo, d(_i) = 4 A) 

- "-salt (+1) counterions {q{+i) 1 eo, 4 A) 

The colloidal surface charge density is = 0.95 eonm^^ = 
0.152 Cm~^, which is in the experimentally relevant regime. 
200 trivalent counterions correspond to a concentration of 
C(+3) — 30 mM. Up to risait = 1300 monovalent salt ion 
pairs were added (concentration Csait = 196 mM). The sys- 
tem is overall electro-neutral and the colloid is fixed in the 
center of the simulation box. 

All particles interact via the Coulomb interaction 

Koulomb('") = 4— (1) 
r 

where the qi are the charges of the respective particle types 
and = 7.1 A is the Bjerrum length in water at room 
temperature. Furthermore, the steric repulsion between 
the particles is modelled via the core interaction 

, II 7" - f off < ?'ci 

L , otherwise 

(2) 

This is the well-known shifted Lennard- Jones potential 
with an offset of Toff. To take into account only the re- 
pulsive part of the potential, we choose rcut = 2^/^ and 
Kihift = £■ Also, we choose cr = 1 A and e = IkT to ensure 
the same "hardness" of all interactions, so that the dis- 
tance of closest approach is well-defined. All particle sizes 
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Fig. 1. Integrated charge Zint against the distance from the 
colloid center r for different concentrations of additional salt 

Csalt . 

are modelled via the values of ros for the different inter- 
actions, which is defined by = ^{di + dj) —a, where di 
and dj are the diameters of the respective particle types. 

The system was simulated using a standard molecular 
dynamics (MD) algorithm using the Velocity Verlet algo- 
rithm in the (N, V, T)-ensemble, with an MD time-step 
of St = 0.01. To speed up the computation of the long- 
ranged electrostatic interactions, the P3M algorithm [27] 
was employed and tuned such, that the maximal error in 
the forces was /Imax — 0.001. The water environment was 
modelled implicitly via the Langevin thermostat (T = 1, 
7 = 0.5). Hydrodynamic interactions were neglected, as 
we are only interested in static equilibrium observables. 

All simulations were performed using the simulation 
package ESPResSo [H]. 

3 Results 

Figure [l] shows a plot of the integrated charge Zint{r) 
against the distance to the center of the colloid for dif- 
ferent concentrations of the additional salt Cgaiti which is 
defined by 

^int(r) = Qo + V / q^p^{r')dr' (3) 
i J" 

where the Pi{r) are the local densities of the different par- 
ticle types in a distance r to the center of the colloid. 
The plot clearly depicts, that in all simulations, the (neg- 
ative) charge Qq of the colloid is overcompensated and has 
a maximum of up to 16eo at a distance of slightly more 
than two ion diameters from the colloidal surface. 

In the following, we define the effective charge Z^s 
of the colloid as the maximum of the integrated charge 
distribution Zi^t- This definition of the effective charge 
corresponds to the static model of Tanaka and Grosberg 

Figure [2] contains a plot of Z^ff and the contributions 
of the different ion types against the concentration Csait 
of additional monovalent salt. Details are revealed by the 
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Additional 1:1 salt concentration c^gi,(mmol/l) 

Fig. 2. EfTective charge Zofi (maximum of the integrated charge 
Zint) against the concentration of additional salt Csait, and con- 
tributions of the different ion types. 
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Fig. 3. Radial distribution function g{r) of the trivalent coun- 
terions on the colloidal surface that contribute to the effective 
charge for different concentrations of additional 1:1 salt Csait- 




plots of the contributions of the different ion types in the 
figure. Keep in mind, that only by subtracting the bare 
charge Qq of the colloid from the contribution of the mul- 
tivalent (+3) ions Zgff (4,3-), the different contributions are 
brought to the same scale. This stresses, that at all values 
of Cgaitj the multivalent ions are responsible for the largest 
part of the compensation of the colloid's bare charge. 

Note, that our definition of ZeS is not identical to 
the definition used in any of the previously cited theories, 
where only the multivalent ions contribute to the effective 
charge. To be able to directly compare with the theories, 
one can look at the plot of the quantity —Qo + Zoff^(+3), 
i.e. the contribution of the multivalent ions Z^ff (_(_3) mi- 
nus the bare charge Qo of the colloid. In contrast to what 
Shklosvkii's and Levin's theories predict, this quantity 
monotonously decreases with growing salt concentration. 
Also, it does not show a maximum, as the theory of Pi- 
anegonda et al would predict [15]. 

However, it is also apparent, that the number of multi- 
valent ions that contribute to Zcs decreases with growing 
concentration of additional salt, up to a point where the 
multivalent ions alone would not suffice to compensate 
Qo anymore. This is in contradiction to a fundamental as- 
sumption of the OCP theories, that we will discuss in the 
following section. 

Instead, the contribution of monovalent (+1) counte- 
rions Zf,s,{+i) rapidly increases and is only partly com- 
pensated by the also increasing (negative) contribution of 
monovalent colons ^cff,(-i)- 

3.1 Strongly correlated liquid but no screening 

All OCP theories predict, that the mutlivalent ions form a 
strongly correlated liquid on the colloidal surface, provided 
that the coupling is large enough. The plasma parameter 
of our systems is I2D — {nasisv'^y^^ — 6.3, therefore we 
are in the strong coupling regime. This view is supported 
by the plot of the radial distribution function g(r) of the 
trivalent ions that contribute to Z^s in figure [3] which 
shows a correlation hole at very low distances and a (weak) 



3 



Csait = Ommol/1 




Csait = 91 mmol/1 Csait = 196mmol/l 

• : (+3)-ions <-: (+l)-ions >-: (-l)-ion 



Fig. 4. Snapshots of the colloid and the ions at the surface 
that contribute to the effective charge Z^g for different con- 
centrations of additional 1:1 salt Csait. 

peak at a distance of about a = 17 A. Also, the simulation 
snapshots in figure |4] show, that the trivalent counterions 
indeed bind and completely cover the colloidal surface and 
form what looks like a liquid- like structure. 

According to Shklosvkii's theory, the only effect of 
adding monovalent salt to the system should be to en- 
hance the Debye-screening of the electrostatic interactions 
between the multivalent ions. This effect should be most 
prominent, when the Debye screening length of the mono- 
valent salt, which is defined by 

rs = (S^^BCsait)"'/' (4) 

is smaller than the preferred distance a between the mul- 
tivalent ions. In our simulations, we used values of the 
screening length of down to rs = S A, which is clearly 
smaller than the preferred distance a between the triva- 
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Fig. 5. Contribution of the different ion types to the charge 
density p against the distance from the colloid center r for dif- 
ferent concentrations of additional salt Cgait • Note the different 
scales on the j/-axis! 



lent ions. Still, g{r) does not show a shift in the position 
of the first peak, which would be expected, if the pre- 
ferred ion distance a would be reduced. This refutes the 
idea, that Debye screening plays a significant role for the 
structure of the multivalent ions close to the surface. 

Instead, we see a small decrease in the height of the 
peak of g(r), meaning that the ions density decreases. This 
observation is consistent with the fact that less trivalent 
ions contribute to the effective charge at higher salt con- 
centrations. 



3.2 Monovalent ions on the surface 

As was shown above, even though the contribution of the 
trivalent ions to the effective charge is most important, 
also the monovalent counter- and colons contribute to the 
effective charge Zeff. Where are those ions located? 

Figure |5] shows a plot of the charge densities qiPi of 
the different ion types against the distance to the colloid 
center r. As was to be expected, with increasing Cgait, the 
density of monovalent ions close to the colloid increases. 
Interestingly, the maximum of the monovalent counterion 
density is at the same distance to the colloid as the max- 
imum of the trivalent counterions for all salt concentra- 
tions. Note, however, that the maximum of the integrated 
charge Zjnt is at much larger values of r « 60 A. 

While the value at the maximum of the density of the 
(+l)-ions increases with growing Cgait, the maximum value 
of the density of the multivalent ions decreases, i.e. some 
of the multivalent counterions are replaced by monovalent 
counterions, which reduces the effective charge. This fact 
shows a weakness in the OCP theories, as these do not 
take into account that monovalent counterions play any 
role for the structure on the colloidal surface. 

On the other hand, the monovalent colons form a weak, 
smeared-out second layer, that can also be seen as an in- 
dicator of the charge reversal of the colloid. 
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Fig. 6. Densities of the different ion types on the colloidal 
surface relative to the corresponding bulk densities Pshoii.i/po.i 
against the concentration of additional 1:1 salt Cgait- 

3.3 Multivalent ions in the bulk 

Figure |6] shows a plot of the quantitity Pshoii,i/po.ii where 
the po,i — Nil are the bulk densities of the different ion 
types and the PshoU,i are the respective densities of the ions 
in the shell that contributes to the effective charge. Not 
surprisingly, at Csait = 0, the density of the trivalent ions 
close to the colloid is much greater than the bulk density 
of the trivalent ions, while the density of the monovalent 
colons is smaller than in the bulk. 

The value of /Oshcii,(+i) / Po.C+i) for the monovalent coun- 
terions is almost constant for increasing salt concentration 
Csait, which indicates that the increasing contribution of 
monovalent (-1-1) counterions to the effective charge is sim- 
ply due to the growing bulk density of these ions. 

On the other hand, the value of Pshcii,i/po,i of the (+3)- 
counterions and the (-l)-coions decreases when more salt 
is added. This shows, that with growing salt concentra- 
tion, more and more of these ions prefer to stay in the 
bulk. 

This can be interpreted in favour of the theory of Pi- 
anegonda et al |15j . who take into account that the (neg- 
ative) solvation energy of the multivalent ions in the bulk 
is larger when more salt is there. 

This effect reduces the reversed charge with growing 
salt concentration, which is also the reason why the theory 
of Pianegonda et al is the only OCP theory that correctly 
describes a decrease in the effective charge at higher salt 
concentrations. 



4 Conclusions 

To conclude, our simulations give support to the basic 
assumption of the OCP theories, namely that the over- 
compensation of the colloidal charge by multivalent ions 
is well described by a strongly correlated liquid. However, 
the theories have a few shortcomings when it comes to the 
description of the effect of additional monovalent salt. 

In particular, the theories predict, that with growing 
salt concentration, the strongly correlated liquid can pack 
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more densly on the colloidal surface, as the multivalent 
ions are screened by the monovalent salt. This assump- 
tion is not confirmed in our simulations. High salt concen- 
trations seem not to influence the size of the Wigner-Seitz 
cell, and consequently, we do not see giant charge reversal. 
The giant overcharge observed in the parameter regime of 
Tanaka and Grosberg [18 is in our point of view due to 
the special choice of parameters. For the strong coupling 
regime considered, all ions form Bjerrum pairs, hence the 
concept of screening becomes meaningless. 

Instead, with growing salt concentration, monovalent 
counterions are pushed onto the colloidal surface, and they 
replace some of the multivalent ions, which lowers the ef- 
fective charge. This effect is not included into the current 
OCP theories. 

Finally, we observed, that not only the trivalent coun- 
terions are pushed out of the region close to the surface 
of the colloid when more salt is added, but also some of 
the monovalent colons leave the region. This is consistent 
with the assumption of Pianegonda et al |15j . that the 
multivalent counterions and the monovalent colons can 
gain favourite solvation energy when more salt is in the 
bulk. 

Therefore we conclude, that the theories of Shklovskii 
and Levin fail to describe the phenomenon of colloidal 
charge reversal under the influence of additional 1:1 salt 
correctly, as they are missing important concepts. On the 
other hand, the more recent theory of Pianegonda et al 
|15j contains the necessary ingredients for a good descrip- 
tion, even though the quantitative agreement with simu- 
lations is still not satisfactory p5] . 
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